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Abstract. 

We show that the non-adiabatic perturbations between Ricci dark energy and 
matter can grow both on superhorizon and subhorizon scales, and these non-adiabatic 
perturbations on subhorizon scales can lead to instability in this dark energy model. 
The rapidly growing non-adiabatic modes on subhorizon scales always occur when the 
equation of state parameter of dark energy starts to drop towards —1 near the end 
of matter era, except that the parameter a of Ricci dark energy equals to 1/2. In 
the case where a — 1/2^ the rapidly growing non-adiabatic modes disappear when 
the perturbations in dark energy and matter are adiabatic initially. However, an 
adiabaticity between dark energy and matter perturbations at early time implies a 
non-adiabaticity between matter and radiation, this can influence the ordinary Sachs- 
Wolfe (OSW) effect. Since the amount of Ricci dark energy is not small during matter 
domination, the integrated Sachs- Wolfe (ISW) effect is greatly modified by density 
perturbations of dark energy, leading to a wrong shape of CMB power spectrum. The 
instability in Ricci dark energy is difficult to be alleviated if the effects of coupling 
between baryon and photon on dark energy perturbations are included. 
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1. Introduction 

The present acceleration of tlie universe, suggested by various observations [H El |3],is 
one of the most important puzzles in cosmology. Possible explanations for this puzzle 
can be based on the modification of theory of gravity or introducing a new form of 
energy whose pressure is negative called dark energy |1] - |H]- Many models of dark 
energy have been proposed in literature from various motivations. The interesting one 
is motivated by the holographic principle [9j. 

The holographic principle implies a relation between the short distance (UV) cut- 
off and the long distance (IR) cut-off in quantum field theory [9]. Thus the amplitude 
of vacuum energy, which may play the role of dark energy, is controled by the long 
distance cut-off whose length scale should be comparable with one of the cosmological 
length scales. Choosing the Hubble radius as an IR cut-off, the present amplitude of 
vacuum energy is in agreement with observations Nevertheless, the vacuum energy 
evolves like matter, so that it cannot play the role of dark energy. The equation of 
state parameter of vacuum energy can be negative if particle or future event horizons 
are used as the IR cut-off [TOl [TT] , but for the choice of particle horizon it cannot be 
negative enough to drive the present acceleration of the universe. So we are left with 
the future event horizon but this choice encounters a causality problem [12]. It has been 
shown that with the inverse square root of the Ricci scalar, which is propertional to the 
curvature radius of spacetime, as the IR cut-off, vacuum energy can be a viable dark 
energy candidate, the Ricci dark energy [I3] where 

P.= f^, (1) 

where pd is the energy density of dark energy, R is the Ricci scalar and a is a constant 
parameter. The reduced Planck mass, Mp = I/a/SttG, is set to be unity. 

In order to study the dynamics of cosmological perturbations for Ricci dark energy, 
we follow the literature [H] to suppose that the above equation is also valid in 
the perturbed spacetime, which should be possible because the relation ([T]) can also 
be obtained from some arguments about quantum gravity effects [13]. In [H], the 
adiabaticity between dark energy and matter perturbations is assumed, i.e. Spm/Spd = 
Pm/Pd, where 6pd, Pm and 6pm is the density perturbations in dark energy, energy density 
of matter and density perturbations in matter respectively. For this case, the instability 
can occur in the future, and can be avoided if a > 1/3. In this work, we will check 
whether the adiabaticity between dark energy and matter is preserved in general, and 
study the dynamics of non-adiabatic perturbations between dark energy and matter as 
well as the instability in Ricci dark energy model. We suppose that the dynamics of 
Ricci dark energy can be presented in terms of the perfect fluid and the conservation of 
energy- momentum. 

In section 2, we derive some useful relations for the background dynamics of Ricci 
dark energy. The dynamics of density perturbations is studied in section 3, and the 
conclusions are given in section 4. 
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2. Background dynamics 

Using the trace of Einstein equation, we can write the energy density and pressure in 
terms of Ricci scalar as 

R = pm + Pd-^Pd- (2) 

Here, pd is the pressure of dark energy and Pm = Pc + Pb, where pc and pb is the energy 
density of Cold Dark Matter (CDM) and baryon respectively. Substituting ([1]) into the 
above equation, we get 
1 



2 

Pm + Pd 1 

a 



(3) 



so that the equation of state parameter of dark energy is given by 



Wd=l(r + 1- -], (4) 



3Hwdr , (5) 



3 \ a 

where r = Pm/ Pd- The evolution equation for r can be obtained from the energy 
conservation of matter and dark energy 

^ _ dr _ Pm f Pm Pd 
dt Pd \Pm Pd 

where if = a/a is the Hubble parameter and a is the cosmic scale factor. Using (jlj), the 
solution of the above differential equation is 

WdAa"^^ 

= : — ^ ■ (6) 

Here, Wd = I — 2/a and 

A-^, (7) 

ro + Wd 

where ro is the present value of r. We will use subscript o to denote the present value 
of the variables through out this paper. Since pm = Pmo'^"^, the evolution of pd can be 
presented in terms of r as 

Pd = ^—a-' . (8) 
r 

In order to obtain the present acceleration of the universe, Wd must be negative, so that 
in the early epoch r is nearly constant and given by 

2 

?^early = -Wd = - - 1 . (9) 

Thus, it follows from (jl]) that dark energy evolves like a matter in the early epoch, i.e., 
Wd ~ 0. 

In figure ([1]), the evolutions of r for various values of a are plotted. From the 
figure we see that when looking back in time r reaches reariy rather quickly, i.e. ([9]) is 
approximately valid during matter domination. Thus, during matter domination, the 
density parameter of matter and dark energy is given by 
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Since a should be in the range 1/3 < a < 1 [14J, (ITU]) shows that Qm is smaller than 
unity and is not small during matter domination. We will see in the next section 
that this feature greatly influences the ISW effect. 
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Figure 1. The evolution of r for various value of a is presented by thick lines, while 
the corresponding is presented by thin lines. The present values of flm and fid are 
chosen such that Omo = 0.26 and ildo = 0.74 respectively. 



3. Evolution of perturbations 



In order to study dynamics of density perturbations, we start by perturbing the trace 
of Einstein equation around FRW background as 



Since we suppose that ([T]) is also valid in the perturbed spacetime, the above equation 
can be written as 



5pd 



5pm + 5pd 1 - 



a 



(12) 



Unlike the background case, the perturbations in CDM and baryon evolve 
differently, so that we have to analyze their effects on the perturbations in dark energy 
separately. We will express all perturbed quantities in Newtonian gauge in which the 
line element is 

ds^ = -(1 + 2^)dt'^ + a\t){l - 2^)6ijdx'dx^ . (13) 

According to the conservation of energy momentum tensor, the evolutions of density 
contrast Sc = Spc/ Pc and velocity perturbations Vc of CDM obey [151 HB] 

S' = 3$' - kVr , 



(14) 
(15) 
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while the evolution equations for density contrast 6i, and velocity perturbation Vb of 
baryon are given by [T71 118] 



61 



3$' - kvb , 

- nvb + A;* + clkSb + Rt~^ (v^ 



Vb 



-Uvb + A;^ + Tc ^ 



(16) 
(17) 



where a prime denotes derivative with respect to conformal time rj = J dt/a, a. subscript 
7 denotes photon, is the sound speed of baryon, R = (4/3)p^/pfe and is the Compton 
mean free path of a photon scattering off a baryon. Here, Tc depends on the baryon 
sound speed and the coupling between baryon and photon. In the early epoch, photon 
and baryon are tightly coupled, so that [H] 



Tc 



V 



1 + R 

Tr 



kcl5b + 



R 



1 + R 



k { -6^ 



a^] +'Hvb + V 



1 



+ k ( - + a: 



R 
1 

4S 



U {vb - 



Tr 



1 + R 
1 + 2'H 



a 
a 

T, 



—Vb + k{ -5^ 



2a. 



(19) 



1 + i?. 

where is the anisotropic perturbations in photon. 

The evolution equations for the perturbations in dark energy can also be obtained 
from the conservation of energy momentum tensor. In the case where the anisotropic 
perturbations in dark energy is absent, one gets [IHl fT6] 



s' = ~m 



+ 3 (1 + Wd) - kud ■ 



n {3wd -l)ud + k— + (1 + Wd)k^ 



(20) 
(21) 



where Ud = (1 + Wd)vd is the momentum density of dark energy. 
Substituting (|T2|) into the above equations yields 

5'^= -Hr (5„ - 5rf) + 3 (1 + Wd) - kud , 



U^ = 'H {?>Wd -l)Ud + k 2 (f^m - 5d) + WdSa 



+ {l + Wd)k<iJ 



(22) 
(23) 



{PcSc + Pb5b)/Pr. 



where 5rr 

known. Thus, we use f|T4|) - (H 

y = 3$' - kv. 



These equations can be solved if the evolution of 5m is 
to derive the evolution equation for 5m as 

(24) 
(25) 



Pb 



Vm= - 'HVm + km + — Tc . 

Pm 

Here, Vm = {pcVc + PbVb)/Pm- 

In our case, it is more interesting to describe the perturbations in dark energy in 
terms of the following variables 

^= {1 + Wd)5m- 5d, U = {1 + Wd)Vm- Ud. (26) 

It can be seen that the adiabaticity between dark energy and matter perturbations 
corresponds to A = {/ = 0. Using (l22i) - (!25i) The evolution equations for A and U can 
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be written as 
A' = nrA - kU . 



U' = n (Swd -l)U + nwd(^^-^^Vm-^(^^-l^A + Wd^ - 4^ (5„ + (1 + Wd) 

From these equations, one can show that on superhorizon scales the evolution of A is 
governed by 

A" +n{l-3wd- r) A' - {n' + n^) rA 

~ - nwd (- - a) kvm - (- - 4^ 5^ - (1 + Wd) — kTc , (29) 

and on subhorizon scales A obeys 

A"+'Hil-3w,-r)A'-'^(--l]A 

~ - nwa (- - a) kvm - wX (- - 4) 5^ - (1 + Wd) — kTc . (30) 
\a J 3 \a J Pm 



Here, we keep the terms on the RHS of ( ]29l) . because in the case of our interest A can 
be very small so that its evolution on large scales can be affected by these terms. 

To study how A evolves, we start with the simplest case where the main 
contributions to the matter perturbations in (fT2!) come from CDM perturbations, i.e., 
Sm ~ and Vm ~ Vc, and the contribution from Tc becomes negligible. In the radiation 
and matter eras, Wd ~ 0, so that the second and fourth terms on the RHS of ( l28l) can 
be neglected. Hence, it follows from fl27|) and f l28|) that if the perturbations in dark 
energy and matter are initially adiabatic in the early epoch, i.e., A = U = 0, they will 
remain adiabatic at later time as long as ^ both on superhorizon and subhorizon 
scales. Nevertheless, if the perturbations in dark energy and matter are not adiabatic 
in the early epoch, the non-adiabatic perturbations can grow both inside and outside 
the horizon. During the radiation and matter eras, the RHS of fl29l) can be neglected, 
so that the solutions of (l29l) are 



^ I CiTj'^ + C2 during radiation domination, 

C3?7^'" + C4?7~^ during matter domination 

It can be seen that A has a growing mode, and the constant adiabatic mode corresponds 
to ci = C2 = C3 = C4 = 0. Similarly, during the radiation and matter eras, (|30l) becomes a 
homogeneous differential equation. To solve this equation, we use the variable x = ^/v^^ 
where p = (r — l)/2 and p = r — 1 during radiation and matter domination respectively. 
In terms of this variable, (!30|) becomes 

k ''"early , ^ f 2 ^ 



x"- 



X^O, (32) 



3 ' T]^ V" 4^ 

where p = rcariy/2 and p = reariy — 1/2 during radiation and matter eras respectively. 
Since krj ^ 1 inside the horizon, the non-adiabatic perturbations grow inside the horizon 
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as 

^ _ f rfi?7('--i)/2eV'"<=--iy/3'='? + 4r^(''-i)/2e-A/''<=-i>-/3fc'? during radiation domination 

I d^'q^'^^ey^^^^^'''^ + c?4r7''"^e" v^^^^'^'"' during matter domination, 

Again, the constant adiabatic mode corresponds to di = d2 = d'^ = d^ = 0, and |A| 
can grow inside the horizon. Since l/3<a<l, risin the range 1 < r < 5 during 
the radiation and matter eras. Thus, from small non-adiabatic perturbations between 
dark energy and matter in the early epoch, |A| and also \U\ can grow by many orders of 
magnitude, and can become unphysically large at late time. This implies an unphysically 
large amplitude of 5d and also 5m together with the metric perturbation due to the back 
reaction which follows from the perturbed einstein equation [16] . 

3H (^(^ + H'^j + ^^= ~ ^ (^'-^^ + + ^'^'^^) 

= - ^ i^rSr + i^m + ^d)5m " ^d^) , (34) 

where a subscript r denotes radiation including photon and neutrino, and a subscript 
m denotes matter including CDM and also baryon in general case. 

Roughly speaking, the existence of the unphysically large density perturbations, 
which implies an instability problem in dark energy, can be avoided if the non-adiabatic 
perturbations between dark energy and matter disappear or have a very small amplitude. 
However, near the end of matter domination, these non-adiabatic perturbations can be 
generated although they disappear in the early epoch. This is because Wd becomes 
negative, so that (!29l) and (!30!) become inhomogeneous. Obviously, the constant 
adiabatic solution A = f/ = cannot be the particular solutions of these equations. 
On large scales, the terms on the RHS of fl29l) cannot influence the evolution of A so 
much because their amplitude is small. Nevertheless, inside the horizon, to create the 
structure in the universe, 5m must grow, consequently the inhomogeneous part of (!30ll 
leads to rapid growth of |A|. As a result, 5d, 5m and the metric perturbation rapidly 
grow and become unphysically large. However, it is easy to see that the inhomogeneous 
part of fl30|) will always vanish if a = 1/2, and hence the generation of non-adiabatic 
mode after matter domination can be avoided for this value of a. These qualitative 
considerations are in agreement with the numerical integration presented in figures ([2]) 
and ([3]) for the case where a = 0.4 and a = 1/2 respectively. In these figure, the initial 
conditions for A and U are chosen such that A = U = 0, and the perturbation variables 
are normalized such that the curvature perturbation ( defined below in ( l39|) equals to 
4. The numerical solution of the evolution equations for the perturbations and also the 
CMB power spectrum are obtained using CMBEASY |T9]. From these figures one can 
see that the instability problem in Ricci dark energy can be avoided if a = 1/2. 

It follows from the above consideration that if the contributions from baryon 
perturbations in (fT2l) can be neglected, the instability problem in Ricci dark energy can 
be avoided if the perturbations in dark energy and matter are adiabatic initially, and 
a = 1/2. Let us now compute the initial conditions for the perturbations in dark energy 
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Figure 2. The evolution of $, (5m and A. The thick hnes represent the perturbations 
mode that enters the horizon about the present epoch, k = 2.4 x lO^'^^Mpc^^, while the 
thin lines represent the perturbations mode that enters the horizon about the time of 
radiation-matter equality, k = 1.6 x 10~^Mpc~^ . For the thin lines, we plot logj^g 1*^1; 
log^o \Sm\ and log^g |A| instead of their actual values. 
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Figure 3. The evolution of $, 6m and A. The thick lines represent the perturbations 
mode that enters the horizon about the present epoch, k — 2.4 x 10^'^Mpc^^, while the 
thin lines represent the perturbations mode that enters the horizon about the time of 
radiation-matter equality, k = 1.6 x 10~^Mpc~^ . For the thin lines, we plot logj^p '^m 
instead of Sm- 



and the other species needed for the numerical integration for density perturbations and 
CMB power spectrum. Since the energy density of Ricci dark energy depends on Ricci 
scalar, the amplitude of the density perturbation in this dark energy depends on the 
metric perturbations. The relation between the density and metric perturbations can 
be obtained from (1T]1 as [Tit [20] 

5pd = -SR = a 2—$ - 4—^ - 3$ - 6H<^> - 9H^ - 6{H + 2H^)<I> , (35) 



2 
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5^ = -___<|.. (37) 



where the vanishing of the anisotropic perturbations in radiation is not assumed. Hence, 
on the superhorizon scales and in the radiation dominated epoch, 

5. = ^ = ^ (h + $ = -^^4?^ - - ^fi. + 2) $ = -2$ . 

Pd Pm+Pd^ ^ \lm + ^ld\ 11) 

To obtain the above equation, we suppose that $ is nearly constant on large scales 
because we use the adiabatic initial conditions for photon and neutrino, and the 
adiabatic perturbations between dark energy and matter A = is the constant mode. 
Using the adiabatic initial conditions derived in [16], the relation between the initial 
value of the density contrast of photon 5-^ and metric perturbations can be written in 
Newtonian gauge as 

10 
5 + 2?]^ 

This implies that in the early epoch if the perturbations in matter and dark energy 
are adiabatic, i.e., 5^ = 5m = —2$, the perturbations in matter and photon cannot be 
adiabatic. This conclusion holds even though the anisotropic perturbations in neutrino 
are neglected, because, in this case, br = —2$ = Sm in the early epoch. 

To study how the non-adiabaticity between radiation and matter perturbations 
influences OSW effect, we derive the evolution equation for the curvature perturbations 
on slices of uniform total energy density by ignoring the anisotropic perturbations 
in photon and neutrino , and treating photon and neutrino as a single perfect fluid. 
Furthermore, since in the early epoch dark energy evolves like matter and 6d = Sm, we 
also treat dark energy and matter as a single fluid whose density contrast^m obeys (^^. 
Here, a subscript m denotes this combined fluid. The evolution of the density contrast 
of the total energy density 6t = Spx/pr, where 6pT = Spr + Sp^ and pr = Pr + Pm, is 
governed by the general evolution equation for the perfect fluid 

S' = -m - Wt] 5t + 3 (1 + Wt) - k{l + wt)vt , (38) 

\6pT J 

where wt = {l/3)pr/pT is the total equation of state parameter, and vt = Vr + Vm is 
the total velocity perturbations. From the definition of the curvature perturbations on 
slices of uniform total energy density [16] 



C=-$ + TT7T^^, (39) 

3(1 + Wt) 

and (|38|) . one can show that on large scales 

C = ^ {5PT - c'Jpt) = ^Pnad , (40) 

Pt+Pt Pt + Pt 

where = (l/3)p^/p^ is the total adiabatic sound speed. It follows from the above 
equation that if the perturbation in radiation and fluid m are adiabatic, Pnad will 
vanish. Consequently, ( is constant on large scales (cee e.g. [21]). In the radiation 
dominated epoch, the main contribution to ( comes from the curvature perturbation on 
the hypersurfaces of uniform radiation energy density, Crd = — $rd + ^rrd/^ = — 3$rd/2, 
where the subscript rd denotes the value in radiation dominated epoch. Similarly, during 
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the matter domination, (^^^ = — $md + (^mmd/S = — 5$md/3, where the subscript md 
denotes the value in the matter dominated epoch. Thus, the constancy of ( yields 
'^'md = 9$rd/10 [22]. In the case of non-adiabatic perturbation, we have to know how 
the entropy perturbations in radiation and rh fluid evolve. According to the conservation 
of energy, 6r obeys 

SI = 4$' - hvr . (41) 

Using the above equation and we get S' = —k{vr — fm), where S = 36r/4: — dm is 
the entropy perturbations between radiation and m fluid. Hence, on large scales S ~ 
constant |22] , so that 5r — 5m until the matter domination. However, the detail analysis 
shows that the non-adiabatic perturbations between radiation and matter can grow in 
early epoch [23]. Since this growth of non-adiabatic perturbations between radiation and 
matter does not influence our rough estimation below, we assume that S ~ constant, 
and therefore 

^ 1 p'r Pt ^ 1 e 1 4 

Pnad = OPr " Opr = -Opr - - Opr , 42 

6 Pj. Pr 6 6A\lr + 6\l^ 

which vanishes in the radiation dominated epoch and then becomes negative when 
Vtr < 1, because in our notations $ is negative and 6pr is positive. This implies that 
( increases during the transition from radiation domination to matter domination, so 
that Crd < Cmd and hence |$md| > 9|$rd|/10, which is in agreement with the numerical 
solution plotted in figure (jl]). Using the approximations Sr ^ Sfn — —2$ and \E' ~ $ on 
large scales during last scattering epoch, one gets 

Go + ^ ~ Oo + $ = ^(5. + $ = ^<I> , (43) 

where Go is the temperature perturbations in Newtonian gauge [22] and the above 
expression is evaluated at last scattering. This equation shows that the non-adiabaticity 
between radiation and matter perturbations enhances the OSW effect compared with 
the adiabatic case in which Qq + = /3. 

The main contributions from the density perturbations in Ricci dark energy to the 
density perturbations in the universe come from the ISW effect which may be represented 
in terms of On large scales and after last scattering, $ is governed by the perturbed 
Einstein equation 

$ = _ ^ ^Q^Sm + nj,) , (44) 

Here, we have supposed that $ ^ \1'. Since $ and the density contrast have opposite 
signs in our cases. The contribution from 6^ can reduce |$| because it cancels the 
contribution from the term — The contribution from 6d increases if Qd during matter 
era increases, and when the contribution from 6^ increases l^l decreases. Usually OSW 
effect is partially canceled by ISW effect [22], a smaller |$| leads to a smaller cancellation 
and therefore a large amplitude of CMB power spectrum at low multipoles. For the case 
of Ricci dark energy, is not small during the matter era, so that the contribution from 
6d can significantly reduce |$|, leading to a large amplitude of the CMB spectrum at 
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Figure 4. The large scale evolution of $, 6d and ftd^d for Ricci dark energy 
model (thick lines) and dark energy model with constant equation of state parameter 
Wd — —0.9 (thin lines). For later model, we use adiabatic initial conditions for all 
species, while forRicci dark energy model we set a = 1/2 and A = U = initially. 



low multipoles. The reduction of |$| due to the contribution from 6d is shown in figure 
dl]). In figure ([5]), the CMB power spectrum for the case where the instabihty problem 
in Ricci dark energy can be avoided, i.e., A = U = initially and a = 1/2, is plotted. 
It can be seen from the figure that the CMB power spectrum is greatly enhanced at 
small multipoles because the SW effect is modified due to the density perturbations in 
Ricci dark energy. 



Co 



12000 
10000 - 

8000 - 

6000 

4000 

2000 




1000 10000 



Figure 5. The CMB power spectrum for the case where a — 1/2 and A = U — 
initially. 



Finally, we include the baryon perturbations in 6pm in ( 1T2|) . so that the term Tc on 
the RHS of ( 128|) cannot be neglected. Using the same consideration as above, one can 
see that this term can generate a growing non-adiabatic mode inside the horizon, leading 
to an instability problem in Ricci dark energy model. Up to our trial, the generation 
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of this non-adiabatic mode cannot be avoided by adjusting the initial conditions for 
the perturbations, special values of a or can be canceled by the other terms, hence the 
instability problem is not easy to avoid in this case. The evolution of Sm and A for 
this case is plotted in figure 
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Figure 6. The evolution of 6m and A. The thick hnes represent the perturbations 
mode that enters the horizon about the present epoch, k — 2.4 x 10~^Mpc~^, while the 
thin lines represent the perturbations mode that enters the horizon about the time of 
radiation-matter equality, k = 1.6 x lO^^Mpc^^ . For the thin lines, we plot log^o |*f |, 
\ogiQ \5m\ and log^Q |A| instead of their actual values. 



4. Conclusions 

We study the evolution of density perturbations in Ricci dark energy model. We 
extend the analysis in [H] by allowing the non-adiabatic perturbations between matter 
and dark energy. In the case where the contributions from baryon perturbations are 
neglected, if the perturbations in matter and dark energy are adiabatic initially, they 
will remain adiabatic both on superhorizon and subhorizon scales until the equation of 
state parameter of dark energy becomes negative. When the equation of state parameter 
of dark energy becomes negative, the non-adiabatic perturbations between matter and 
dark energy perturbations will be generated even though they disappear in the early 
epoch. This non-adiabatic mode is a growing mode that can lead to an instability in 
Ricci dark energy, i.e., the density perturbations in dark energy, matter, and metric 
perturbations grow extremely large. The generation of this non-adiabatic mode can 
be avoided if a = 1/2 and the perturbations in dark energy and matter are adiabatic 
initially. Nevertheless, the perturbations in matter and photon cannot be adiabatic 
in the early epoch, so that the OSW effect is modified. The main influence of the 
perturbations in dark energy on the density perturbations in the universe is the great 
modification of the ISW effect due to a large amount of dark energy during matter era. 
This modification of ISW effect greatly enhances the CMB spectrum at low multipoles. 
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When the contributions from baryon perturbations are included, the couphng between 
baryon and photon can generate growing non-adiabatic mode that leads to the instability 
in Ricci dark energy. In this case, the generation of growing non-adiabatic mode cannot 
be avoided by simply adjusting the initial conditions for the density perturbation or 
spacial choices of a. 
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